Persistence in Advection of Passive Scalar 



OO 

o 
o 

(N 
-t— > 
O 

O 

(N 

Si 

o 

B 

i 

> 

•4-* 



o 
o 



> 

(N 

in 
o 



OO 

o 



X 



D.Chakraborty * 
Department of Theoretical Physics, 
Indian Association for the Cultivation of Science, 
Jadavpur, Kolkata 700 032, India. 

We consider the persistence phenomenon in advectecd passive scalar equation in 1-dimension. The 
velocity field is random with the (v(k, u>)v(— k, — id)) ~ \k\~ < ~ 2+a K In presence of the non-linearity 
the complete Green's function becomes G" 1 = — iid + Dk 2 + E. We determine E self-consistently 
from the correlation function which gives E ~ k , with /3 = (1 - a) j2. The effect of the non-linear 
term in the equation in the 0{t 2 ) is to replace the diffusion term due to molecular viscosity by an 
effective term of the form Eofc' 3 . The stationary correlator for this system is [Sech(T/2)] 1/ '' 3 . Using 
the self-consistent theory we have determined the relation between (3 and a. Finally, IIA is used to 
determine the persistent exponent. 



The problem of persistence or survival [l[ has attracted 
a lot of interest in the last decade. The persistence prob- 
ability has been obtained both analytically and numeri- 
cally for a large class of stochastic process, Markovian as 
well as Non-Markovian. The random walk problem, dif- 
fusion problem, surface growth, Ising model with Glauber 
dynamics are only few to name (H] |24j. The persistence 
probability or the zero crossing probability is simply the 
probability that the local field <fi(x,t) has not changed 
sign up to time t. For single particle systems such as the 
random walker, which are also Markovian in nature, the 
persistence probability is easy to calculate since the sta- 
tionary correlator of such a process decays exponentially. 
For many body systems where the field cj> has a space 
dependence the calculation of the zero crossing probabil- 
ity becomes complicated. The problem is now two fold- 
first we have to write down an effective equation for a 
single site process by solving the underlying dynamics 
of the many-particle system and then from this effective 
single site equation we have to find out the persistence 
probability. Even though the first part is achieved, the 
second part of obtaining the persistence probability is 
notoriously tough since the resulting single site process 
becomes non-Markovian. 

The simplest of such a process which one can think of 
is the diffusion equation dt<fi = Dd 2 <fi. The fact that this 
is effective single site equation can be seen from the so- 
lution <f>(x,t) = J dx'G(x — x', t)(j){x' , 0), where G(x,t) is 
the Green's function for diffusion equation. The prob- 
lem of persistence in a diffusion equation has already 
been addressed by Majumdar et.al [2j|. They consid- 
ered the diffusion equation with random initial condition 
4>(x, 0) taken from a Gaussian distribution. The two time 
corrleation function C(t\,t2)oi the normalized variable 
X = <j){x, t) I \J (<j) 2 (x, t)) takes the form 



C(t 1 ,t 2 ) = (X(i 1 )X(i 2 )) 



[4<iW(ti+< 2 ) 2 ] D/4 , (I) 



where D is the dimension of space. Now if we make 
the transformation T = \xit, the correlator C(Ti,T2) be- 
comes f(\Tt - T 2 \), with f (T) = [sech(T/2)] D / 2 , which 
is clearly stationary. The stationary correlator for the 
effective single site process is not exponentially decaying 
and therefore the calculation of the persistence exponent 
becomes difficult. The fact that the correlator is not ex- 
ponentially decaying indicates that the effective single 
site process is non-Markovian because of the interaction 
with nearest neghibour sites. Given this stationary non- 
Markovian correlator it then remains to determine the 
persistence probability. Two methods have been devel- 
oped to address this problem, the Independent Interval 
Approximation (IIA) Q and the "series expansion" Q 
approach. In this present article we will use IIA to eval- 
uate the exponents. 

Knowing the information about the persistence expo- 
nents for a diffusive process, it is natural to ask what 
would be zero crossing probability when the diffusive pro- 
cess is augmented by an advection term. For a simple 
diffusive process, if L be the relevant length scale (say 
the size of the container), then the time to diffuse to a 
distance L is simply t\ = L 2 /2D. If, however, the par- 
ticles are advected then the the time for them to diffuse 
through a distance L is r 2 = L/v, where v is the advec- 
tion velocity. The ratio of the two time scales is 



n/^ = 2S/Re, 



(2) 



where 5* is the Schmidt number and Re is the Reynolds 
number. The Schmidt number is of the order of unity and 
therefore, it follows that the mixing time due to advection 
is smaller than the pure diffusive process. We, therefore 
expect that the exponents will be greater than those for 
the pure diffusive process. 

The advected passive scalar equation reads 



at 



+ v(x, t) ■ = DW 2 



(3) 



together with 
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V -v = 0. 



(4) 
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In 1-dimension however, the constraint imposed by 
Eq.(j4]) is relaxed. Instead, in 1-dimension, we will con- 
sider a random velocity field drawn from a given distri- 
bution. In particular, the velocity- velocity correlation is 
given by 

(v(x, t)v{x',t')) = 2Dg{\x - x'\)8{t - t'). (5) 
The Fourier transform of Eq. ([3]) in one dimension is 

^ +i J " P ' t] = Dk2 ^ [kl t] (6) 

and the velocity-velocity correlation transforms to 

(v(k, t)v(k', t)> = 2Dg(\k\)S(k + k')S(t - t'), (7) 

where (f>{k,t), v(k,t) and g(k) are the Fourier transform 
o{(j>(x,t), v(x,t) and g(x) respectively. We choose <?(|fc|) 
as a power law decaying function, that is, 



9{\k\) = 



|fc|(a+a)- 



(8) 



with < a < 1. In absence of the non-linear term the 
Green's function for Eq.© is 



GW 1 = -iw + Dk 2 



(9) 



The effect of the nonlinearity is to replace the zeroth 
order Greens function by its most general form 



G- 1 = -iu; + Dk 2 + Z = G n 



Hence, 



G 



Go 



1 + G E 



G (l - EG + E 2 G 2 



(10) 



(11) 



The correlation function ((p(k, oj)(p(—k, —u>)) can be writ- 
ten as 

(4>(k, cj)4>(-k, -u)) ~ GG = G Go-Go£Go+G £ 2 G§+... 

(12) 

We will determine the self-energy E self consistently from 
the correlation function. 

Taking a Fourier transformation of Eq. © in time do- 
main we get 

- iuj(f>(k,uj) + i J ^^~Pv( k ~ ~ u')<i>(jp,w') 

= -Dk 2 4>(k,uj) (13) 
while the velocity corrleation function becomes 
(v(k, u)v(k', </)) = 2Dg(\k\)S(k + k')6(uJ + uj r ). (14) 
We now make a perturbative expansion in <fi and write 



= <p Q + e<pi + e (p 2 



(15) 



Substituting this in Eq. lfTB"]) . the zeroth order solution is 
4> {k,Lu)[-iuj + Dk 2 ] = Mk,0) 

Mk, w) = ; = G o(fc, u)Mk, 0) (16) 

[ — ilo + L>k z \ 

In the first order the solution for <f>x{k,ui) is 

(f>i(k,w) = G (k,w)[-i / 2~^T pv(k-p,u}-u)')<j> (p,u)')] 

(17) 

pv(k-p,w- w')0i(p,w')] 



while the solution for cf>2 becomes 

dp du>' 



4>2{k,u) = G (k,cj)[-i 
= G (fc, lo 



2tt 2tt 
dp du>' 



' I -^^-pv(k - p,uj - uj')G (p,uj') 

dqduj" , „ - , 

-I I —-^qv(p-q,uj -u )(f>o{q,UJ ) 



2tt 2tt 



(18) 



To evaluate S self consistently we need to cal- 
culate the corrleation function {4>\{k, u})<j)\{— k, — u>)) 
(<f>2{k, lu)4>o(— k, — lu)). We assume that the non-linear 
contribution to the total Green's function G will domi- 
nate over the Dk 2 term [25]. Hence, we rewrite Eq. fTU)) 
as 



G- 1 



-iu + E, 



(19) 



which shows that ui and E have the same dimension. 
The correlation (4>2(k,uj)(f)o{— k, —to)} is then 



i{k,u))(j> Q (-k, -u)) = -G (k,(j){ 



dp da/ 



-i I Y~-^-Qv(p- q,u' - u")(t>o{q,u") \<j> Q (-k,-u>) 



pv{k -p,u>- w')G (p, uj') 

> 

(20) 



A little algebra simplifies the above expression to 



i(k,cj)4>o(—k,—u>)) = —Go(k,u) 



dp du>' duj' 1 



2ir 2n 2tt 

kp(v(k — p,u> — uj')v(p — k,u> — uj"))Go(k, w")Gq(— k, — u>) 



(21) 



The velocity- velocity correlation gives a 5(u> — u>") which, 
after the oj" integral becomes 



b 2 (k,u))<f> (-k, -co)) = -G (k,u>) 

Go(p,w') 



dp dui' kp 
2^27|fc-p|( 2 +") 

G (fc,w)Go(-*,-w) 

(22) 
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We now turn our attention to (<pi(k, a;)^>i(— k, — u>)) 
which is given by 



The correlation C(t\,t2) = (<p(x,ti)(f>(x,t2)) for a fixed x 
is given by 



i(k,uj)4>i(-k, -u>)) = -G (k,uj) 



dp duj' dq duj' 1 



2tt 2tt 2tt 2tt 
pq(v(k — p, uj — uj')v(—k — q, —uj — uj')) 



C(h,t 2 ) 



(4o(p,u/)$o(q,u>")) 



G (-k,-u) 

(23) 



dfc(<KMi)<K-M 2 )) 



I[£ (t 1+ i 2 )]-V/3 



(31) 



The velocity-velocity correlation introduces a 5{p + 
q)5(u) / + oj") while the average ((j>o(p,uj')(j)o(q,uj")) gives 
us S(p + q). Integrating over the q and uj" variable we 
get 



Define the normalized variable X(t) = 
cj>(x,t)/-y/{(t> 2 (x,t)). Then, the correlation {X(ti)X(t 2 )) 
in terms of C(ti, t 2 ) becomes, 



C(i 1; i 2 )= (X^Xfe)) = C(ti,«2)/VC(tl,ti)C(t2,t2) 



1 (fc, w)0i(-fc, -w)) = G (fe,w) 



dp duj' 



P 



(*i+*a) 



2tt 2tt |fc-p|( 2+Q ) 



G (p,a;')Go(-p, -w') 



G (-*,-w) 



Making the usual transformation lni = T, Eq. 
comes 



(32) 
be- 



(24) 



The second term in Eq. (fT2|) has the same structure of 
(<^>i(fc, ui)(f)i(— k, — u>)) while the third term has the same 
structure as (4> 2 (k, u))4>q(— k, — ui)). Thus Eq. (f2"2"|) gives 
us 



C(T U T 2 ) 



; 1/2(T!-T 2 ) + e -l/2(Ti-T 2 ) 



1/0 



= [SechC^-^)] 17 ^/^-^!) 



dp do/ kp 
2^ 2tt |fc-p|( 2+Q ) 



G (pV) (25) 



E ~ /fed"")/ 2 
while from Eq. (f2"4")) we get 

dp dtu' p 2 1 



2tt 2tt \k - p\( 2 + a ) uj' 2 + D 2 p 4 



(26) 



(27) 



(33) 

The correlator in Eq. (|33p is now stationary since it de- 
pends only on the difference T 2 \ and non-Mar kovian. 

To determine the persistence exponent we adapt the 
method of IIA as explained in Ref Q. The basic as- 
sumption is that the intervals between the successive 
zeros of X(T) are statistically independent. We will 
briefly outline the method here. The first step is to 
construct the variable a — sign(X). The correlator 
A(T) = (o-(T)ct(O)) is given by 



Since uj ~ S, and neglecting Dp 2 term compared to E, 
power counting yields 



A{T) = - arcsin[/(T)] 



(34) 



fc(l-a) 



(28) 



which gives us the same result as in Eq. (|26[) . We remark, 
in passing, that the result obtained in Eq. (|2l)|) can also be 
obtained by introducing noise term in Eq.([6]). It should 
be noted that for a Kolmogorov like velocity field, a — 
-1/3. 

Thus, the effect of the non-linearity in 0(e 2 ) is to re- 
place the term Dk 2 by an effective diffusion term that 
looks like EqAt. We can, therefore, rewrite Eq.© as 



If p n (T) be the probability that an interval of size T 
contains n zeros of X(T), P(T) be the distribution of 
intervals and Q(T) be the probability that the left and 
right of the interval contains no zeros, then 

Pn (T) = (T)- 1 f dTidT 2 dT„Q(Ti)P(T 2 - Ti).. 

Jo 



.P{T n - T n _i)Q(T - T n ) 



dj)_ 
~dt 



together with 



(29) 



with /? = (1 — a)/2. The two time correlation function 

((p(k,uj)(j)(—k,—uj)) becomes 



(35) 



(36) 



n=0 



$(Mi)£(-fc,ta)> 



_ -£ fc' 3 (tl+t2) 



(30) 



where (T) = —2/A'(0). Taking a Laplace transform of 
Eq.® and using the fact that P(s) = 1 - sQ(s), P(s) 
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and Q(s) being the Laplace transform of P(T) and Q(T) 
respectively, we arrive at 

P-( g ) = -^[1- P^rP^H*) n>l 
(T)s 2 

= 7^[(T)-1 + P( S )]. (37) 
(T)s 2 

Finally, substituting Eq. (|3"T)) in Eq. (f3"rJ)) and carrying 
out the summation over n yields P(s) in terms of A(s), 
that is, 

P(s) = {2-F(s)]/F(s). (38) 
where F(s) is given by 

F(s) = l + ^[l-sA(s)}, (39) 

A(s) being the Laplace transform of A(T). For large T, 
Po(T) ~ e~ eT means that the exponent 9 is given by the 
pole of P(s) or the zero of F(s). 

In our present case, the correlator A{T) is given by 

A(T) = - arcsin[Sech(T/2)] 1/ ' 3 . (40) 

7T 

which gives (T) = iry/4/3. The function F(s) has the 
form, 

2 f°° 

F(s) = 1+ttJ]3s[1 s / dTe- sT arcsin([Sech(T/2)] 1/ ' 3 )] 

f io 

The zeros of the function F(s) are determined nu- 
merically. As a check for numerical verification we took 
the values a = —3,-1 and — 1/3, which corresponds 
to 1/(3 = 1/2,1 and 3/2 respectively. These values of 
1/(3 correspond to the normal diffusion in D = 1,2 and 
3 respectively. The exponents determined numerically 
using these three values of 1/(3 agrees well with the 



exponents reported in Ref[]. Finally, we have taken 
various values of a between and 1 and have obtained 
the roots of F(s) numerically. The obtained values of 
the exponents are listed below in Table A. 



TABLE A 



a 


(3=(l-a)/2 


e 


0.1 


2.22 


-0.29341041 


0.2 


2.50 


-0.312802995 


0.3 


2.86 


-0.336107784 


0.4 


3.33 


-0.364881293 


0.5 


4.00 


-0.401726555 


0.6 


5.00 


-0.451442543 


0.7 


6.67 


-0.524308324 


0.8 


10.00 


-0.64860239 



In conclusion, we have considered the persistence phe- 
nomenon in advectecd passive scalar equation. In 1- 
dimension the velocity is drawn from a random distri- 
bution with (v(k, u)v(-k, -lj)) ~ |fc|-( 2+Q ). The effect 
of the non-linearity is to replace the the zeroth order 
Green's function by it's general form G _1 = —iuj + Dk 2 + 
E, with £ ~ k 13 . We have determined the scaling form of 
E using self-consistent theory, which gives (3 — (1 — a)/2. 
Thus, in C(e 2 ), the non-linearity replaces the original 
dynamics with an effective equation where the diffusion 
term due to molecular viscosity by a term of the form 
Sofc' 3 - We have calculated the two time correlation for 
the effective process which has the form [S~\ J ((T/2)] 1 ^ . 
Finally, we have used IIA to calculate the persistence 
exponents. 
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